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We present an innovative method to resum infrared and collinear logarithms appearing in 
distributions of jet observables in QCD. The method, based on a general master formula with 
applicability conditions, allows resummations at next-to-leading logarithmic accuracy in an 
automated way. As a sample application we present resummed results in hadronic dijet events. 



1 Jet observables in QCD 

A special aspect of QCD in the Standard Model is that the regime where the theory is strongly 
coupled is within reach of current experiments, and indeed unavoidable. This gives rise to 
a number of non trivial phenomena, for instance one encounters non-convergent perturbative 
(PT) expansions and large non-perturbative (NP) corrections even at scales which are, at least 
formally, within the domain of perturbation theory. One of the main aims of today's QCD 
studies is the comprehension of all these phenomena. Apart from its fundamental interest, a 
detailed QCD description is vital for any precision/search physics in modern colliders. 

The study of jet observables, event-shapes and jet-rates, is particularly informative in that 
these quantities are calculable in perturbation theory with a remarkably high accuracy, but 
they are still sensitive to NP, low energy physics. For brevity, in the following I will refer to 
event-shapes only; these are infrared and collinear safe observables describing the topology of 
the final state in high energy collisions. 

In recent years their mean values and distributions have been used to measure the cou- 
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pling constant the colour factors 1 ^ and to study their hadronisation corrections*^ 1 The most 
discriminatory studies make use of distributions. Integrated distributions £(i>) are defined by 
requiring that the value of the observable V(k\, . . . , k n ) - a function of all n momenta - be less 
than a fixed value v 
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These distributions are characterised by two kinematical scales, the hard scale Q of the process 
(for instance the c.o.m. energy in e + e~ collisions, the virtuality of the gauge boson in DIS 
processes, etc.) and an observable scale vQ whose magnitude is determined by the typical 
momenta of secondary emissions. In the more exclusive kinematical region v <C 1 there is 
a strong ordering between these two scales, so that the distribution Eq. ^ contains (large) 
logarithms of the ratio of the two scales L = ln(l/i>). They compensate the smallness of a s , and 
to make predictions one needs to resum logarithmically enhanced (next-to-leading NLL) terms 
to all orders, while in the region where v = O (1) fixed order (next-to-leading NLO) predictions 
are reliable. 

It is interesting to compare today's effort required to obtain NLO and NLL predictions. 
To compute distributions at a given PT order one needs to know the emission probabilities 
(the squared matrix elements) at that perturbative order. Given the QCD Lagrangian, these 
quantities are calculable in principle at any order, though technically it becomes extremely 
difficult to go beyond NLO. Since PT matrix elements are independent of the observable, they 
can be implemented in a general way in Fixed Order Monte Carlos^ These are then widely 
exploited for experimental studies, the only additional input needed being the definition of 
the specific observable in the form of a computer routine. Obtaining NLL descriptions also 
requires the knowledge of the emission probabilities. These quantities mix different logarithmic 
orders, so that one needs to reorganise the PT series so as to account for {resum) all leading 
(LL) exp{a™L™ +1 } and next-to-leading logarithmic (NLL) exp{a"L n } terms. Also, since one 
is resumming terms to all PT orders, one needs to know the value of the observable given an 
arbitrary number of secondary emissions. These are the reasons why, so far, the description 
of each observable has required a separate analytical calculation^ this limits the experimental 
study of resummed predictions, as compared to fixed order ones. 

2 Automated resummation 

The aim of this project is to automate the NLL resummation of jet observables, so as to ease 
their experimental analysis. To achieve this goal one has to understand the origin of all NLL 
enhanced terms in observable distributions. This reduces to two steps: one needs to know the 
behaviour of the observable when only one soft-collinear (SC) gluon is present in the final state 
(single emission properties); furthermore one needs to understand the way in which all emissions 
coherently determine the value of the observable (multiple emission properties). 

2.1 Single emission properties 

To describe the behaviour of the observable due to a single emission we fix an (arbitrary) Born 
event and simulate the emission of one soft gluon collinear to each hard parton (leg) t. We 
assume that the observable can be parameterised as 



where kt, rj denote respectively the transverse momentum and the rapidity with respect to the 
emitting hard leg and is the azimuthal angle with respect to a plane. Our numerical code 
verifies automatically the assumed parameterisation and determines the coefficients ai,bi, and 
the functions <?£ (</>). This knowledge allows us to compute all the NLL terms which account for 
hard collinear emissions, soft large angle radiation and inclusive gluon splittings (running of 
the coupling). This gives rise to a single (or simple) distribution S s (u) = e~ R "( v \ where R s (v) 
is a LL Sudakov exponent (complete formulae can be found elsewhere)P' It also happens, by 
construction, that S s is the full NLL resummed distribution for a simpler observable V s whose 
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Table 1: Single emission properties for T m ,A, 1 and 2 label the incoming legs, 3 and 4 the outgoing ones. 
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value is determined only by the largest emission 

V s (h, ...,k n ) = max[y(fci), . . . , F(fc n )] . (3) 

2.2 Multiple emission properties 

To fully resum the observable at NLL one also needs to relate the value of the observable to 
the momenta of all secondary emissions, i.e. to understand the observable specific mismatch 
between the full observable V and its simplified variant V s . An analytical treatment of this 
requires a detailed insight in the kinematics of the recoiling hard partons and, for jet-rates, of 
the recombination procedure, but this effect can be computed numerically in a general way^ 
and is given as a pure NLL function !F(R'), with R' = dR s (v s ) / dln(l/ v s ) . 

2. 3 The master formula and its applicability conditions 

The knowledge of the single and multiple emission properties is combined in a master formula^ 

= Z s (v s )F(R') . (4) 

The formula Eq. |1] does of course not apply to all final state observables in QCD and is accom- 
panied by a list of conditions the observable has to fulfill it should be recursively infrared and 
collinear safef" should vanish in the Born limit; should be continuously global^ should behave 
as assumed in Eq. [5] when only one SC emission is present; additionally it should, at the present 
stage, satisfy some minor (technical) requirements. While this might seem a long list practically 
the limiting condition is the requirement of globalness, while all other conditions are satisfied by 
all observables resummed so far. Most importantly an essential feature of our code is the ability 
to verify all properties automatically and to resum the observable only when the correctness of 
the result is guaranteed at NLL accuracy. 

2.4 Sample output: first ever resummations in hadronic dijet production 

The most interesting application of our method is the resummation of observables in hadronic 
dijet events. As an example we present here results for the indirectly global thrust minor" 
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here pu denote the momenta transverse to the beam, p x i the momenta out of the plane defined 
by the beam-axis and the transverse thrust axis and A = O (1) denotes a rapidity cut. The 
program verifies all applicability conditions; tabulates the single emission properties, see Tab.^ 
computes for all colour configurations the multiple emission function J 7 , see Fig. ^ (left); and 
exploits the master formula Eq. |1] to obtain the full NLL resummed prediction. Fig^ (right) 

"Note that the second term guarantees the continuous globalness of the observable. 
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Figure 1: The function T (left) and the NLL resummed distributions (right) for the indirectly global thrust 



shows differential distributions, for several underlying hard subprocesses, for the Tevatron run 
II regime, ^fs = 1.96 GeV. We select events having two outgoing jets with Et > 50GeV and 
\r)\ < 1, fix A = 1, use the CTEQ6M parton density set^ corresponding to a s (Mz) = 0.118 and 
set the factorisation and renormalisation scales \xp and /ir at the partonic c.o.m. energy. 

2.5 Final remarks, conclusions and outlook 

Results presented here have the quality of analytic NLL predictions, though no human in- 
tervention is needed. This implies that any hadronisation model can be applied, studies of 
renormalisation and factorisation scale dependence can be carried out and, since the answer is 
free of spurious subleading terms, a matching with fixed order is feasible. 

The first step in this project has been the numerical computation of multiple-emission effects, 
which allowed us to resum three observables in e + e~-collisions at a time^ We have now a code 
which fully automates the resummation of a large class of jet observables; as a result we obtained 
the first ever resummations in hadronic dijet production. As a next step we aim to automate the 
matching with fixed order results. This will open up the possibility to carry out a vast amount 
of phenomenological studies. 

Acknowledgements 

I am enjoying working at this project with Andrea Banfi and Gavin Salam. A special thank to 
Yuri Dokshitzer and Pino Marchesini for stimulating conversations and precious suggestions. 

References 

1. S. Bethke, J. Phys. G 26, R27 (2000). 

2. S. Kluth et al., Eur. Phys. J. C 21, 199 (2001) and references therein. 

3. P. A. Movilla Fernandez, S. Bethke, O. Biebel and S. Kluth, Eur. Phys. J. C 22, 1 (2001). 

4. D. Graudenz, |arXiv:hep-ph/9710244| J. Campbell and R. K. Ellis, Phys. Rev. D 65 (2002) 
113007 ; Z. Nagy, Phys. Rev. Lett. 88, 122003 (2002). 

5. A. Banfi, G. P. Salam and G. Zanderighi, arXiv:hep-ph/0304148 and in preparation. 

6. A. Banfi, G. P. Salam and G. Zanderighi, JHEP 0201, 018 (2002). 

7. M. Dasgupta and G. P. Salam, Phys. Lett. B 512, 323 (2001); JHEP 0208, 032 (2002) 
and references therein. 

8. J. Pumplin et al. JHEP 0207, 012 (2002). 



